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i-^ ' Abstract 

f^ . This is a brief review of some of the basic concepts of perturba- 

D I tive QCD, including infrared safety and factorization, relating them 



to more familiar ideas from quantum mechanics and relativity. It is 
intended to offer perspective on methods and terms whose use is com- 
monplace, but whose physical origins are sometimes obscure. 
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1 Asymptotic freedom in QCD 

We begin with a short portrait of quantum chromodynaniics, the unbroken, 
nonabehan gauge theory SU(3). The classical Lagrange density of QCD can 
be represented schematically by pLj 

where q labels quark fields, of mass r/ig, and where G^jy,a[^] is the nonabelian 
field strength, including the self-couplings of the gluon field. We can think 
of this expression as an analogy to electrodynamics, the sum of kinetic terms 
for the quarks and gluons, supplemented by various local interactions. QCD 
is the Yang-Mills gauge theory [2] of quarks and gluons, in which gluons 
are like photons with charge, so that the gluon field is a source for itself. 
This nonlinearity, of course, is part of what makes QCD, and the strong 
interactions it describes, the source of such varied phenomena. It was realized 
early on that the quarks of QCD provide just the right currents to couple to 
electromagnetic and weak interactions, so that previous results based on the 
analysis of those currents ( "current algebra" ) could be taken over essentially 
unchanged [3]. In addition, this theory has just the right kind of forces: the 
QCD charge is "antishielded" , growing larger with increasing distances over 
which it is measured. This is its famous property of asymptotic freedom [1] . 

Let us sketch how the asymptotic freedom of QCD is established. Working 
conceptually, imagine that we definethe strong coupling, g{h/T), as just the 
amplitude for a quark to emit a gluon within a sphere of radius cT, with c 
the speed of light. So defined, g{h/T) is a running coupling, varying with 
the scale at which it is measured. The concept is illustrated in Fig. [H In a 
sense, we send a quark into the sphere, wait a time of order T, and see if it 
comes out accompanied by precisely one gluon. The amplitude for this to 
happen is given by an infinite set of perturbative diagrams, each representing 
a particular set of quantum mechanical histories. We show some of the lowest 
order diagrams in Fig. [H 

Now the diagrams within the sphere are described by integrals that do 
not converge, because there are simply too many states with one or more 
additional gluons at very large energy. Nevertheless, with a bit of work, we 
can compute the T-dependence of g{h/T), or in more familiar notation, the 
/i = /i/T-dependence of as(/i) = 5'^(/i)/47r. With Uf different (flavors of) 



g(h/T) = 




Figure 1: The running coupling defined by a sphere of radius cT . 
quarks, we find 

where we adopt the common notation 60 = 11 — 2nj/3, and where a(/io) 
and the scale Aqcd can be thought of as integration constants. The value 
of Aqcd is set by any boundary condition for as(/io) at any scale /^o- In 
other words, it is set by nature, once we learn to measure as at a given scale. 
Eq. ([2]) expresses asymptotic freedom, according to which as(/i — *• cxd) — *• 0. 
In QCD, the colors of virtual gluons "line up", like neighboring magnets, a 
feature that depends on both the spin and the self-interactions of the gluons 
built into the Lagrangian. The smaller the sphere, the fewer the lined-up 
magnets, and the weaker the interaction. 

An essential result of asymptotic freedom in QCD is that radiation be- 
comes weaker as momentum scales increase, or equivalently distances (like 
cT above) decrease. In effect, near a color source, the coupling constant is 
weak, a feature that leads to the famous approximate "scaling" observed in 
deep-inelastic scattering, which we will describe below. Correspondingly, far 
from a source, the coupling constant appears to grow, a feature that is at least 
consistent with (although by no means ensuring) the observed confinement 
of colored quarks and gluons. 

In the years leading up to the discovery of QCD, a template O El [7] had 
been developed to connect the behavior of the running coupling in any field 
theory with what we now call parton distributions, which we will denote 



as /i///(^,/u), for partons i carrying momentuni fraction ^ of hadron H [S]. 
Here, /z is a renormalization scale, very much analogous to h/T above, and 
can be thought of as determining the scale at which we probe hadron H to 
count these partons i, leading to probability density fi/H{^,fJ')- This probe- 
scale dependence is encoded in sets of "anomalous dimensions", which can 
be computed as power series in the couplings. 
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77v(as) = — 7}^^^ + --- , (3) 

where for QCD, we know that as(/i) vanishes as /i increases. For moments of 
the parton distributions, f{N) = Jq d^^^~^f{0^ ^^e general analysis gives 

mun\ 



fi/H{N,jj) = fi/HiN,no) exp 



1 f^ dfJL 



/2 



2 J, I /i'2 
and with ^^(yu) = 47r/6oln(/i^/AQCj3), we get: 



7jv(as(/^')) 



St^'V^O 



(4) 



Once the 7Ar's were computed at one loop (and eventually all the way to 
three loops [9J), it all worked. 

To get back to the parton distributions, fi/ni^, f^), we can invert the mo- 
ment transform. We can then compute structure function F2{x,Q'^), which 
describes deeply inelastic scattering in terms of the variable x = Q^/2p ■ q 
and the momentum transfer Q^ = —q^ > (see Eq. ( IT6l) below). The data of 
Fig. [2] shows exactly a pattern predicted by the explicit forms of the 7Ar's: ap- 
proximate scaling (Q-independence) at moderate x and pronounced evolution 
(Q-dependence) for small x. Perfect scaling for the structure functions would 
follow for vanishing coupling. This corresponds to //-independent parton dis- 
tributions, as in the parton model, which provided a successful description of 
the first moderate-x data, in which Q-dependence seemed weak if not absent 
altogether [8]. 

As has been widely recognized, the asymptoic freedom of the QCD run- 
ning coupling is a result of historic significance. This is as much because it 
opens the door to new studies, as because it explained previously mysterious 
features of nature. A tongue-in-cheek analogy that I like is 
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Figure 2: (a) Approximate scaling at moderate x. (b) Pronounced evolution 
for smaller x. Data from the HI experiment at HERA. 



In its explanation of approximate scaling (and the violation, or "breaking" 
of scaling), asymptotic freedom is a beginning, not an end. For Newtonian 
gravity, the immediate challenge to the inverse-square law was the three-body 
problem (moon-sun-earth, for example). For QCD it is how to study a theory 
in which the fundamental degrees of freedom are masked by confinement. The 
ultimate goal might be expressed in a similar spirit as 

Nuclear Physics Chemistry 

QCD " QED ■ ^ ^ 

A short summary of questions we must ask in this context include: can we 

• Study the particles that give rise to electroweak currents (quarks)? 

• Study the particles that provide the forces (gluons)? 

• Expand in the number of gluons (i.e., use perturbation theory)? 

In QCD the fundamental quanta are confined, and (at least in the absence 
of extreme temperature and pressure) observed hadrons are bound states. 
The scattering of bound states confronts us with the complex structure of 
these hadrons, and the strong forces that hold them together on length scales 
comparable to 1/Aqcd- A question that was raised often in the early days 
of QCD was quite simply, "Does this make sense at all?" 



2 Learning to Calculate with the Theory 

Our first observation is that not all is hopeless. Certain quantities even in a 
confining theory are quite "perturbation theory-friendly" . 

2.1 Correlations and the S-matrix 

The classic examples of quantities closely related to perturbation theory 
are correlation functions between color-singlet currents at short distances, 
schematically, 

{N\J{z)J{0)\N) = CN{zfi,as{fi)) 

= C^(l,a,(l/z)) , (8) 

in some state N. When N is the vacuum, the primary example is the total 
e"'"e~ annihilation cross section, for which J is the electroweak current. In 
this case, the function Cn can be expressed as a power series in the coupling 
evaluated at the momentum scale of inverse distance 1/z (here treated as a 
simple scalar). Any such quantity, which depends only on one or more short 
distance scale, is said to be infrared safe. When \N) is a nucleon state, such 
matrix elements are related to deep-inelastic scattering. In this case, the 
function Cn is somewhat more complex, and the matrix element is not itself 
infrared safe, but its dependence on the short length scale is still computable, 
using the factorization formalism that we review below. 

Calculating an 5'-matrix element in perturbative QCD, however, is pretty 
hopeless, 

{B ont\A in) = / [Q/fi, m/fi, ^^(/i)) 
= /(l,m/Q,a,(g)) 
= /(Q/m,l,a,(m)), (9) 

where A and B are hadronic states, Q is a hard scale, and m denotes various 
soft scales in the theory, including the masses of light quarks, the (perturba- 
tively vanishing) mass of the gluon, and the strong-coupling scale Aqcd from 
Eq. ([2]). No matter what choice we take for the scale in the running coupling, 
we encounter large ratios of the energy to fixed mass scales. If S-matrix el- 
ements are not accessible, are we doomed to compute only correlations of 
currents? The answer turns out to be "not quite", and here we can turn to 
another strand of the story. 



2.2 Structure of final states: Cosmic rays to quark 
pairs 

As it turns out, not being able to compute S-matrix elements is not the same 
as begin forbidden to "look inside the final state" . In fact, as we now know, it 
is possible to see in certain final states a direct portrait of quarks and gluons 
in the form of "jets" of nearly collinear high-energy particles. The story of 
the term "jets" actually begins before QCD, in fact even before the paper of 
Yang and Mills. This is the tale of particle jets in cosmic rays. 

While tracing back some references, I was surprised to read in a paper 
from 1957, that "The average transverse momentum resulting from our mea- 
surements is pt=0-5 BeV/c for pions [a table] gives a summary of jet events 
observed to date ..." [10]. Evidently, the jets associated with perturbative 
QCD did not by themselves give rise to the term. What was being reported 
in this paper was a spray of particles of very high energy but limited trans- 
verse momentum (a BeV is a GeV), observed in cosmic ray events, as seen 
in emulsions. 

Somewhat over ten years after Ref . [TO] , accelerators had been developed 
to study the multi-GeV annihilation of electrons and positrons into virtual 
photons, which can then decay into anything that carries charge. In the 
meantime, the quark model had been invented, and the quarks carry charge. 
So, what was going to happen? If (as everyone suspected) we wouldn't see the 
quarks because of confinement, what would we see? Jets? In Physical Review, 
Drell, Levy and Yan [TT] took the step of extending the parton model from 
deep-inelastic scattering to e+e~ annihilation, and built into their model the 
same limited transverse momentum that had been observed a decade earlier 
in cosmic rays, describing the limitation as a cutoff: "Because of our cutoff 
kmax "^ 1 1? I • • • The distribution of secondaries in the colliding ring frame will 
look like two jets ..." 

Now this was a real prediction for the nature of final states in e"'"e~ — >■ 
hadrons, and following the spirit of the parton model, Drell, Levy and Yan 
suggested that the angular distribution of the jets would follow the same 
angular distribution as a quark pair, or any other spin-1/2 pair, in Born 
approximation, 1 + cos^ 6, with 6 the angle to the beam axis (see Fig. [3]). In 
this picture, partons "fragment" into hadrons. Whether this would happen 
was a question to ask of both nature and of QCD. Would the final states 
look like this? 

In nature, they did, as shown by the analysis of Hanson et al. at SLAC in 
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Figure 3: Representation of the conjecture that there is a relation between 
quark pairs and jets for e"'"e~ annihilation. 

1975 [I2]. And, in the fullness of time, that's what happens in deep-inelastic 
scattering, in e+e~ annihilation and in hadron-hadron scattering. Figures 
show nature's answer to the question of whether jets exist. 
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Figure 4: A jet in deep-inelastic scattering. Event recorded at the HI exper- 
iment at HERA. 



2.3 How to calculate jet cross sections 

Clearly, we can observe the jets, but we still have to ask whether we can calcu- 
late anything about them. Here we can hope to benefit from the asymptotic 
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Figure 5: A jet pair in e+e annihilation. Event recorded at the Opal exper- 
iment at LEP. 

freedom of QCD, but as we've seen, we have to be careful - the S-matrix 
cannot be treated by short- distance analysis alone. If so, how can we hope 
to compute cross sections? 

We can get insight into the challenges involved, and their possible solu- 
tions, by recalling the related "infrared problem" of QED, and its "solution" . 
As is often the case, the problem is related to asking the wrong question, and 
the solution to identifying the right one. 

In QED, typical exclusive cross sections have infrared divergent correc- 
tions in perturbation theory, which show up as logarithmic dependence on 
the (vanishing) mass of the photon. This happens as soon as we go to the 
order oem = e^/47r correction of a Born cross section, say in electron-electron 
scattering at momentum transfer Q: 

Cr^eLee (Q, rXle, ITL^ -* 0, Oem) ~ "EM ^ee^ee(QMe) In , (10) 

with Pee^eeiQ/'nT'e) a fuuctiou that is finite for vanishing photon mass, m^. 
Following the famous Bloch-Nordsieck analysis [T3], however, we trace this 
divergence to asking an unphysical question, the probability for one or more 
charged particles to scatter, and in the process be accelerated, while emitting 
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Figure 6: Jet pair at a hadron collider. Event recorded by the DO experiment 
at Fermilab. 



no radiation at all. In effect, we are computing the probability of something 
that never happens. 

The classical theory demands radiation, and classical radiation requires 
an essentially unlimited number of very low-energy photons. This is Bohr's 
correspondence principle between quantum and classical mechanics. Rather 
than count the number of photons (zero, one . . . ), Bloch and Nordsieck [13] 
counseled that we introduce an energy resolution, eQ with e ^ 1, and then 
sum over final states with arbitrary photon emission, as long as the total 
energy comes in below the energy resolution. Experimentally, this is not a 
choice, but a necessity, because our apparatus will always miss some photons 
if they are soft enough. At first, however, it sounds complicated. How can 
we sum over all soft photons? But this will not be necessary. 

Following the Bloch-Nordsieck procedure, suddenly the full order oem 
correction with an energy resolution becomes finite by itself, the log of m^/Q 
being replaced by a log of e, as the result of a cancellation between the final 
states with and without an extra photon. As long as Q;EMln(l/e) <ti 1 (which 
is almost inevitable), the correction is small. 



=(1) 

'ee^ee+X(e) 



(Q, me, eQ, Oem) ~ OEM /3ee^ee(Q/"^e 
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The magic (and beauty) of this is that we don't have to sum over an infinite 
number of soft photons, even though this is the root cause of the problem! 

Now let's think about QED in the very high-energy limit. On a closer 
look, we find that the function /3ee^ee(Q/'"^e) itself has a log of Q/me- Given 
that the dominant term depends only on the ratio, we can as well trade 
the high-energy limit for the double, photon-and-electron, zero mass limit. 
But in this case, our energy resolution is not enough to produce finite cross 
sections. If, however, we can solve this problem in QED, we may be able to 
solve it as well in high-energy QCD, where the high-energy limit also involves 
logarithmic enhancements in ratios of momenta to all particle masses. 

We wish to look for quantities that are capable of measurement, have a 
single hard momentum scale, Q, and which nevertheless have no powers of 
ln((5/m), only at worst (m/Q) ln((5/m). Such quantities become functions 
of only that single hard scale, and are calculable as a power series in the 
coupling as(/x), with fi = Q, without introducing any large ratios when 
Q ^ oo with fixed masses, or equivalently m ^ at fixed Q. These are 
quantities for which asymptotic freedom can be naturally applied, and are, 
in the terminology mentioned above, infrared safe. 

We've already seen that an energy resolution alone is not enough for in- 
frared safety. Progress can be made, however, by an analogy to the argument 
for an energy resolution based on the unobservability of arbitrarily soft pho- 
tons. We can just as well say that two exactly coUinear massless particles 
cannot be distinguished from a single massless particle of the same total mo- 
mentum (and total quantum numbers), whether that momentum is soft or 
not. That is, if p^ = and p'"^ = and if p^ and p''^ are collinear, then 
{p + p'Y = as well. So whether the combination is a single particle or two 
particles is not easy to distinguish. 

This approach works, and enables us to take the zero-mass limit for all 
particles in QED and QCD. Roughly speaking, any cross section with an 
energy and an angular resolution is infrared safe in e"'"e~ annihilation [T3] . 
If two particles are closer together in direction than some angle, 6, then we 
treat them the same way as we do a single particle. We'll also see that 
hadron-hadron cross sections of this sort, while not themselves infrared safe, 
contain an infrared safe factor that we can isolate. 

The conditions for infrared safety may also be rephrased in a more general 
form as follows. Any cross section that sums over all states that (1) differ by 
the emission or aborption of soft particles, or (2) by the splitting or recom- 
bination of exactly collinear particles, is infrared safe [HI [T51 [ISl ^7\ [TSl [12] 
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(or contains an infrared safe factor.) It is worth noting that to prove infrared 
safety for jet cross sections requires an extension of the beautiful theorems 
that apply to fully inclusive transition probabilities ^U\. This involves a 
careful reanalysis of perturbation theory, and is especially dependent on how 
the gauge invariance of the theory manifests itself [16] . 

The most direct application of infrared safety is to jet cross sections in 
e"'"e~ annihilation, exactly of the type illustrated schematically by Fig. [3] and 
in experiment by Fig. [51 So long as the resolutions are large, we can represent 
the use of asymptotic freedom for jet cross sections by an appropriate choice 
of renormalization scale, taken here as the total energy Q, 

a {Q/fi, e, 6, ^^(/i)) = a (1, e, 6, as{Q)) , (12) 

where e and 6 represent the energy and angular resolutions mentioned above, 
or more generally other parameters that define the cross section. Computed 
in this fashion, there is no need for a transverse momentum cutoff of the sort 
invoked in Ref . [TT] . The infrared safety of the observable ensures that high- 
er radiation is suppressed by factors of asipr)- Such radiation is present, 
of course, but it influences the infrared safe quantity through calculable cor- 
rections, just as the effects of soft gluons influence QED cross sections in a 
finite way at higher orders. 

2.4 The field-theoretic content of infrared safety 

Summarizing, we recount the "sorrows" of QCD perturbation theory, and 
how they can be overcome, at least in part. First, there is color confinement, 
which may be interpreted as the statement that matrix elements like 

y'd^xe-*P-^(0|T[g,(x)...]|0), (13) 

in which we take the Fourier transform of a quark or other field with a 
nontrivial color representation has no p^ = m^ pole in a Green function, with 
T time-ordering. (This is confinement.) Second, poles at physical particle 
masses, such as p^ = m^ for pions, 

fd^xe-'P-''{0\T[n{x)...]\0), (14) 

are not accessible to perturbation theory. 
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Despite all this, we are able to use infrared safety and asymptotic free- 
dom for such quantities as the total cross section for e"'"e~ annihilation into 
hadrons. What are we really calculating? Totally inclusive examples like 
these are related by the optical theorem to forward-scattering amplitudes of 
the general form 

/'d^a;e-''?-"(0|T[J(x)J(0)] |0) , (15) 

involving color singlet currents, just as in Eq. ([8]). Deep- inelastic scattering 
involves hadronic matrix elements rather than the ground state fTfi [T8] , 

W^, = 2lm^ fd^xe-''>-^<N{p)\T[J^,{x)UO)]\N{p)> 
on J 

+ f Pm - ^M^ j f P- - ^-^ j ^ ^2(x, Q') , (16) 

for electroweak currents, J^, and nucleon states, \N{p)), with x = p ■ q/Q^ 
and g^ = —Q^ < 0. The Fi{x, Q"^) are the same structure functions as shown 
in Fig. [2J For such matrix elements, we will apply factorization properties, 
which will enable us to isolate the infrared safe factors referred to above. 

Another class of color singlet matrix elements enables us to describe jet- 
related cross sections ETl. These look like 



lim i?2 f d^^ fdhf{n)e-''^-y{0\ J{0)T[niTo,{xo,Rn)J{y)] |0) , (17) 

R^oo J J 

with Toj the energy momentum tensor, and n a vector on the unit sphere. 
Such a matrix element represents the action of a calorimetric detector, which 
measures energy flow, and matrix elements such as these are what we really 
calculate when we compute jet cross sections. For a general cross section, 
we introduce a "weight", given by function f{n). As long as all the deriva- 
tives, (F f{n)/dff, of the weight are bounded, individual final states may 
have infrared divergences, but they cancel in the sum over collinear split- 
ting/merging and soft parton emission, because these transitions respect en- 
ergy flow |T6]. We regularize the divergences dimensionally (typically) and 
calculate the long-distance enhancements in amplitudes, only to cancel them 
in infrared safe cross sections. It is this intermediate step that makes many 
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calculations tough, and is part (not all) of why higher-order calculations are 
so difficult. It may be worth noting that one of the goals of a collider exper- 
iment is remarkably similar - to control late stage interactions of particles 
once they enter the detectors. 

3 Extracting Infrared Safety: Factorization 

Any cross section with one or two hadrons in the initial state has inescapable 
long-distance behavior, because a semi-inclusive sum over initial states is 
simply not a practical option. In effect, we can choose the energy of the 
nucleon(s) that initiate our scattering process, and sometimes their spin, but 
little else. By construction, then, cross sections at hadronic colliders are not 
infrared safe. The technqiue of factorization, however, enables us to isolate 
and extract infrared safe dependence in a large set of otherwise long-distance 
phenomena. Here we review the physical basis of factorization, and show how 
the factorization of a process also leads to useful information on its energy- 
dependence, including the evolution of the moments of parton distributions, 
as in Eq. (jl]) above. 

3.1 Factorization 

The general form of a factorized cross section (here multiplied by Q^ to make 
it dimensionless) , is 



QVphys(Q, m) = cusD(Q//i, a.(/i)) ® /ld(;U, m) + O (l/Q^) , (18) 

where as shown on the left, the "physical" cross section a depends generically 
on a hard scale Q and on a wealth of soft scales, denoted collectively by m. 
The soft scales include in general the gluon mass, which is zero, as well 
as various quark masses, and the scale of the perturbative coupling, Aqcd, 
encoded in the running coupling, Eq. ([2]). 

On the right of Eq. flTSl) . we give the schematic factorized form of a, in 
which the Q-dependence and m-dependence are separated. There is a short- 
distance function usb, which is infrared safe, and a long-distance function 
/ld; which for hadronic initial states is not calculable in perturbation theory. 
The short- and long-distance functions are linked by a convolution, denoted 
0. For deep-inelastic or hadron-hadron scattering the convolution is in par- 
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tonic momentum fractions, "x" , which is transformed into a simple product 
by the moments leading to Eq. (jlj) above. 

Dimensional analysis requires that we introduce a new scale, fi, the fac- 
torization scale, so that usb and /ld can be nontrivial functions of Q and m, 
respectively. In effect, the factorization scale marks the boundary between 
short- distance and long-distance dependence. 

As indicated in Eq. (ITS]) , factorization is not normally an exact result, 
but it often holds up to corrections that behave as inverse powers of Q. For 
many important examples, such as unpolarized deep-inelastic scattering cross 
sections, corrections enter only as 1/Q^, and are negligible for many purposes 
once Q reaches several GeV. 

In the most familiar examples, including the structure functions Fi and F2 
in flT6|) . the /ld are parton distributions, and we shall refer to them as such. 
The parton distributions themselves can be expressed in terms of expectation 
values [23] in hadronic states that fix light-cone components of the momenta 
of the partons in question. We take the light-cone components for any vector 
v'^ as f^ = (l/v^)(t'° ± v^), with v"^ = 2v^v~ — v^. A vector whose only 
nonvanishing component is v^ or v~ is light-like, v"^ = 0. 

For example, the (spin-averaged) distribution of quark q in nucleon A^ 
with momentum p^ = p'^6^^, and spin s, is 

fM^,t^') = JeT V-"""''"'" <^(P'^)Ulo^i/">o±) 

z ^ J-00 Ztt 

X ^^+<!>^{y~,0)q{0)\N{p,s)> . (19) 

We can compare this form to the matrix element for currents, Eq. fITBl) . In 
this case, the factorization scale, /i, enters because we must renormalize the 
product of quark fields that are separated by a light-like distance y~ in the 
minus direction. The operator $(?/", 0) is a "gauge link", between the two 
fields, whose purpose is to render the matrix element gauge invariant, and 
which is defined by 



<!>n{y ,0) =Pexp 



ry 
-ig / dl n- AUn'') 
Jo 



(20) 



with n^ = 6fj,-. Here the field A^ = J2a ^a'^a is given as a matrix in terms of 
the relevant generators of SU(3), which for quarks are in fundamental (3 x 3) 
representation. The symbol P denotes an ordering of these color matrices 
along the path between / = and I = y~. 
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3.2 Prom factorization to evolution 

If we can factorize a cross section as in Eq. flTSl) . its Q-dependence is calcula- 
ble. As such, we can compute it systematically in extensions of the standard 
model that include new heavy states, which modify the short-distance behav- 
ior of the theory. "New physics" , then, is embedded in a calculable fashion 
in cusD. While not calculable in perturbation theory, the functions /ld are 
"universal" , portable from one factorizable process to another. 

The key to the portability of parton distributions, is their "evolution" , 
which enables us to compute their dependence on the factorization scale 
\25\ . Calculable evolution is not a separate assumption, but rather a direct 
consequence of the factorization in Eq. (1181) . We need only observe that the 
physical cross section cannot depend on the factorization scale, 

= /i— lncrphys(Q,m). (21) 

We can thus separate dependence on Q and m by requiring that the /i- 
dependence of the short-distance function cancel that of the long-distance 
function, 

din USB p, , ^^ ^Iii/ld .r>o\ 

^^-- = -P{a.{^)) = -/^^^ . (22) 

The "separation constant" P{as) can depend only on those variables that 
the short- and long-distance functions hold in common: the coupling and 
the convolution variables. Eq. fl22l) is an evolution equation. We can solve 
it to relate parton distributions at one yU to another, and therefore, because 
we can always choose /x = Q in Eq. flTSl) . we can relate the cross section 
at one Q to that at another scale, up to corrections associated with the 
expansion of usBiQ/l^, cUsIa*)) in the strong coupling. Of course, this analysis 
requires that a^ remain small in the range over which we wish to evolve. 
Schematically, then, we can exhibit the cross section's dependence on the 
momentum transfer as 

crphys(Q, m) = ujsbC^, c^siQ)) ® exp <^ - J —^P (as(/w')) \ /ld(Qo, m) , 

(23) 
just as in Eq. (j3]) for the moments of structure functions. 
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3.3 The pattern of a factorized cross section 

A large class of hadronic cross sections can be factorized, as long as they 
are defined in a manner consistent with the energy flow interpretation de- 
scribed in the previous section. This involves, in general, observing a jet-like 
structure in the final state and summing over soft radiation between the jets. 
The general structure of any such observable falls into a form that can be 
represented schematically as 



dcra+b^jets{i}iQ) _ r r jja'+b'^{di} ci^'+b'^{di} JT J 

-TT;^ - J a' /a ® h'/b ® ^IK ^ ^KI ^ H •Jd, ■ 

"V jets{i} 

(24) 



We can think of this expression as recounting a (quantum-mechanical) story: 
evolved incoming partons represented by fa' /a, fv/b collide and exchange mo- 
menta at short distances. A function Hjk describes quantum corrections at 
that scale (Q), where K and / identify color exchange in the amplitudes and 
their complex conjugates, respectively. These indices are in a color tensor 
basis that reflects the numbers and color representations of all the "active" 
partons, a', 6', and {di\ [26]. In general, the color exchange at short dis- 
tances influences the development of the system at long distances, through 
a color-exchange-dependent soft function. Ski describing the production of 
soft particles. Finally, the production of energetic particles and jets is de- 
scribed by a set of functions J^^, each specifying the fragmentation of the 
parent parton di of jet i. These fragmentation processess are mutually inco- 
herent, with a universal evolution into the final states that is itself the result 
of this factorization. 

Eq. flMl) holds in general to all powers of the coupling, with power cor- 
rections in hard scales. The latter, however, can be quite complex, involving 
ratios of the maximum soft energy to jet energies: Esoft/Ejet, but also inverse 
powers of the energy of soft radiation. That is, we also anticipate "power cor- 
rrections" of the form m/Esoft, with m any of the long-distance mass scales 
in the theory [27] . On the perturbative level, the very presence of a factoriza- 
tion involving soft, jet and short- distance functions ensures more elaborate 
evolutions, involving double-logarithmic corrections [28] . 

It is worth noting that the original cosmic ray jets were not of this sort. 
Their properties are not computable in quite the same way, because for the 
most part they lack truly high-momentum transfer subprocesses, represented 
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by Hki in Eq. (12^ . For recent applications of perturbative QCD to such 
"inclusive" proton- nucleus and nucleus-nucleus cross sections, see Ref. |25] . 

A generalization of Eq. flMl) applies in hadronic scattering to high-transverse 
momentum (pr) single-particle inclusive cross sections. In this case, the jet 
functions of Eq. fl24l) are replaced by fragmentation functions, 

doA+B^H+x{PT) = Yl daA+B^c+xipr/zfi) (g) Dh/c{z, rUc, fJ,) 

c 

^0{m\lvV) , (25) 

with a sum over fragmenting partons c. The cross section da includes parton 
distributions for the initial state. Here, following the formalism developed by 
Collins and Soper [2^, the fragmentation function, Dh/c can be defined as 
a vacuum expectation value similar to those above for the distributions, but 
now involving creation and annihilation operators for the observed hadron. 
For a gluon to fragment to hadron H, for example, the function is [21] 

with G^u the gluon field strength and a'jj- the creation operator for hadron 
H. The relevant ordered exponential, or gauge link, for this process is 

<|.(f'^-')(x-) = Pexp 



oo 



(27) 



where n^ is a lightlike vector in the opposite direction to the jet. For gluon 
fragmentation, the gauge field A^"''^^^ is an 8 x 8 matrix in the adjoint repre- 
sentations of SU(3) generators. Such a gauge link gives rise to a nice set of 
diagrammatic rules, in terms of "eikonal lines" in the x~ direction {n^ = 5^-), 
with vertices —ign'^ x [group factors) and (linear) propagators i/{n ■ k), il- 
lustrated in Fig. [71 To the jet, as it fragments, all that's left of the rest of 
the world is a gluon source moving in the opposite direction, whose entire 
influence is summarized by the eikonal line. Similar considerations apply to 
the parton distributions, Eq. (TT9l) . 

3.4 The classical basis of factorization 



Where do factorized cross sections like Eq. fl2^ come from? In the following 
we review an argument based on the classical Lorentz transformation proper- 
ties of fields, and point out the subtleties of gauge fields in particular [30] . An 
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1/(n-k)^ 



Figure 7: Graphical representation of the ordered exponential, Eq. ([2 

argument based on a classical picture isn't far-fetched, precisely because, as 
noted above, the correspondence principle is the key to the origin of infrared 
divergences. Any accelerated charge must produce classical radiation, and 
infinite numbers of soft gluons are required to make a classical field. Thus 
the classical field has a lot to tell us about the radiation of soft partons. 

Having said this much, we consider the situation illustrated in Fig. [8], in 
which one "bound state", approaches another (from the left in the figure) 
at a relativistic velocity /3 — >■ 1 in the Xa-direction, carrying with it various 
point charges, its "partons" . The coordinates of the bound state on the left 
are indicated by unprimed variables, those on the right by primed. We will 
refer to the former as the projectile, the latter as the target. 

Suppose the partons of our projectile are sources for a massless scalar 
field, whose magnitude we denote by (p. In their own rest frames, the sources 
produce a simple l/\x\ potential. By definition, the magnitude of a scalar 
field at any point in space-time is independent of the coordinate system in 
which it is observed. We can thus start with an expression for our massless 
scalar field in the rest frame of the projectile, and simply reexpress it in terms 
of the coordinates of the target. To do so we use x^ = 'j{Pct' — x'^) = 7A', 
where as usual 7 = (1 — /3^)^^/^. This gives 

where g is a charge and xt the distance of closest approach, which is trans- 
verse to the motion. Naturally, the field is maximized in the target coordi- 
nates at the time of closest approach, where A' = 0, that is, at t' = 4^x'^. 
At this target time, the magnitude of the field is simply q/x'rp. At all other 
values of the time t', however, the field of the oncoming projectile partons 
is proportional to an explicit factor of I/7. In summary, the scalar field 




■,y,z,t 



Figure 8: Schematic representation of the field of an oncoming particle. 

transforms "like a ruler", that is, at any fixed A' 7^ 0, the field decreases 
like 1/7 = Vl — (3'^ as 7 — *■ 00. This is to say that for any fixed time in 
the target frame before closest approach, the field of the projectile decreases 
rapidly as the velocity of the projectile approaches the speed of light. This 
is just a consequence of length contraction in elementary special relativity. 
When an observer riding on the projectile (!) measures a distance X3, then 
an observer sitting on the target measures a much larger distance. 

Next, we suppose that the sources of the projectile couple to the electro- 
magnetic field instead of a scalar field, producing in their own rest frames the 
same l/|a;| potential, but now as the zeroth component of the vector A^(x). 
The following array compares gauge fields to scalar fields from this point of 
view. We compare, on the one hand, the A^ component of the field in the 
projectile frame to the A^ component in the target frame, found by Lorentz 
transformation, and on the other hand the longitudinal (third) component 
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of the electric field in both frames, 

field X frame x' frame 



scalar 



-3- 9 

|x| (x2,+72A2)i/2 



gauge A\x) = ^ A'%x') = ^^.^^%y,. 

field strength E,{x) = ^ E'^{x') = ^,.^;:^2%-sf. ■ 

We can ask the same question of the electromagnetic potential and field 
strength that we posed for the scalar field: at a fixed target time before the 
point of closest approach, how does the field observed at the target depend 
on the velocity of the projectile? The answers for the gauge field and the field 
strength are strikingly different. The gauge potential is actually independent 
of 7 as /5 — > 1 for any fixed A' ^ 0! The vector potential (at least its time 
component) is not contracted at all. On the other hand, the field strength, 
as represented by E^, decreases as 1/7^, which is a much more rapid decrease 
than even the scalar field. 

These two behaviors are, of course, consistent, and are reconciled by the 
realization that the vector field of a relativistic charge approaches a total 
derivative in the target (primed) frame as /3 — ;> 1, 

A' ^{x') = g-^ In {(3ct' - 4) + 0(1 - /3) . (29) 

ox^ 

The bulk of this A'^ is actually an unphysical polarization, and can be re- 

— * 

moved by a gauge transformation. In contrast, the physical "force" field E 
of the projectile does not overlap the target until the moment of closest ap- 
proach. "Advanced" effects in the electric field are corrections to the total 
derivative in Eq. fl29l) . and hence are of the size 



^ 2 



1 r / irl^ m'^ 



y^^ 



2uj 



/2 



(30) 



where m is the mass of the projectile, and u' its energy in the target frame. 
This is a power-suppressed behavior, and a typical initial-state correction to 
factorization. 
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Factorization expresses this contraction effect. As the oncoming projec- 
tile approaches the speed of hght, the appearance of its field is essentially 
instantaneous at the time of closest approach. The projectile then cannot 
affect the internal structure of the target, or vice- versa, and the target's in- 
ternal structure is thus effectively universal among all projectiles, so long as 
the latter are sufficiently relativistic. The initial state structure of the target 
and projectile can then both be summarized by multiplicative factors, and 
these are the parton distributions of Eq. (^^. 

This argument, of course, applies only to initial-state interactions, signals 
exchanged before the hard interaction. For final-state processes to respect 
factorizations like Eq. (1241) . it is also necessary that we define the observable 
in a manner consistent with infrared safety. In addition, for factorization to 
hold, we must require that there be a hard scattering. Otherwise there is no 
well-defined time at which the scattering occurs, and indeed no sharp dis- 
tinction between the initial state and the final state. If there is a well-defined 
hard scattering, however, low-momentum transfers after that scattering are 
too late to affect the large momentum transfer process(es), such as the cre- 
ation of jets or of heavy particles. Similarly, the fragmentation of partons 
into jets of hadrons is too late to know details of the hard scattering, leading 
to the factorization of fragmentation functions. 

3.5 Factorization in perturbation theory 



Perturbative arguments for factorization in QCD [221 [23] are, unfortunately, 
much more complex than the simple classical pictures above. Nevertheless, 
the physical observations we have just made have a direct correlation in 
perturbation theory, which is worth pointing out. We consider a soft gluon, 
of momentum k emitted by a fast quark, whose momentum p'^ is on-shell 
{p^ = rn^) just after this interaction. In perturbation theory, this will be 
associated with a factor like 



u{p) {-rg^n (//fe)/,'!^ = -(^) (-^^) ^ + (IR f^^^t-) 



(31) 



where in the second form we have used the Dirac equation, and have sup- 
pressed terms proportional to k, which are infrared finite, as well as color fac- 
tors. In an arbitrary perturbative diagram, the vector p'^ on the right-hand 
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side will be contracted with the propagator of the soft gluon that carries 
momentum k. Now suppose we were to choose a gauge for which p ■ A = 0, 
in which case the gluon propagator is given (with p^ = 0) by 

In this gauge, the soft gluons decouple from the quark. This argument can 
easily be generalized beyond lowest order, and applies to the entire set of 
collinear partons, whether quark, antiquark or gluon, in that jet. No gauge 
choice like this, of course, can decouple soft gluons from more than one jet 
at a time. But the existence of such a gauge for each jet implies that soft 
gluon couplings cannot resolve more than the direction and overall color of a 
jet [231 [31] . This is the origin of the "universality" of soft gluon interactions, 
and their summary in terms of eikonal lines like those of Eq. ( 1271) and Fig. 
[71 which play a central role in factorization for perturbative QCD. 

4 Conclusion 

We have summarized a few of the major results of perturbative QCD, which 
underly the basic applications of the theory to hadron-hadron and hadron- 
lepton collisions at large momentum transfer. We have presented justifica- 
tions wherever they can be found, in both classical and quantum intuition. 

The coming decade will see unprecedented applications of the ideas and 
methods of perturbative QCD at the Large Hadron Collider, in proton- 
proton, proton-nucleus and nucleus-nucleus experiments. Whether as a pesky 
background to new physics searches, or as a subject of interest in its own 
right, QCD, with its self-generated scales and evolving degrees of freedom, 
will remain a benchmark for our understanding of physics at its most chal- 
lenging. 
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